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Abstract 

New emerging distributed multimedia applications provide guaranteed end-to-end quality of 
service (QoS) and have stringent constraints on delay, delay-jitter, cost, etc. The task of QoS 
routing is to find a route in the network which has sufficient resources to satisfy the constraints. 
Multi-constrained routing is difficult because different constraints can conflict with one another. 
In particular, the delay-cost-constrained routing problem is NP-complete. 

We propose a heuristic algorithm for this problem. The idea is to first reduce the NP- 
complete problem to a simpler one which can be solved in polynomial time, and then solve the 
new problem by either an extended Dijkstra's algorithm or an extended Bellman- Ford algorithm. 
When the extended Dijkstra's algorithm is used, the total time complexity of the heuristic 
algorithm is 0(x 2 V 2 ); when the extended Bellman-Ford algorithm is used, the total complexity 
is 0(xVE), where x is an integer defined solely by the algorithm. We prove the correctness of 
our algorithm by showing that a solution to the simpler problem must also be a solution to the 
original problem. The performance of the algorithm can be adjusted by tuning the value of x. 
A larger value for x results in a higher probability of finding a solution and a higher overhead. 

1 Introduction 

Quality of Service (QoS) routing has been attracting considerable attention in the research commu- 
nity recently [5, 9, 10, 11, 12]. The routing requests are typically specified in terms of constraints. 
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For example, a bandwidth constraint requires the bandwidth of a routing path to be not less than a 
given lower bound, and a delay (cost) constraint requires the total delay (cost) of a path to be not 
greater than a given upper bound. The multi-constrained routing problem is difficult because dif- 
ferent constraints can conflict with one another. In particular, the delay-cost-constrained routing, 
i.e., finding a route between two nodes in the network with both end-to-end delay and end-to-end 
cost bounded, can be formalized as a multi-constrained path problem (MCP), which is NP-complete 
[4, 12]. 

Many existing QoS routing algorithms [9, 10, 11] consider a single constraint. The competitive 
routing strategy discussed in [9] considers only the bandwidth constraint. The Salama's algorithm 
[10] and the Sun's algorithm [11] consider the delay constraint, while using heuristic approaches to 
minimize the cost of the found route. Some multi-constrained routing algorithms work for special 
cases where the constraints are not independent. For example, in the Ma's algorithm [5], delay, 
delay-jitter and buffer space are all functions of bandwidth, which makes the problem solvable in 
polynomial time. Jaffe [4] proposed a distributed algorithm to solve the NP-complete MCP problem 
in pseudo-polynomial time. A more in-depth discussion of the related work can be found in Section 
3. 

We propose a heuristic algorithm for the MCP problem with a polynomial time complexity. 
The idea is to first reduce the NP-complete problem to a simpler problem which can be solved in 
polynomial time by an extended Dijkstra's or Bellman-Ford algorithm, and then solve the simpler 
problem to find a solution path. We prove that a solution path for the simpler problem is also 
a solution path for the original problem. The performance of the algorithm is predictable and 
adjustable. It is predictable in the sense that when certain condition is satisfied the algorithm is 
guaranteed to find a solution. It is adjustable in the sense that the probability of finding a solution 
can be increased when a special integer value defined by the algorithm is increased. 

The rest of the paper is organized as follows. In Section 2, the heuristic algorithm with two 
constraints is presented first; it is then generalized for an arbitrary number of constraints; the 
delay-cost-constrained routing algorithm is discussed next; finally, the performance of the routing 
algorithm is studied by experiments. The related work is covered in Section 3. Section 4 draws the 
conclusion. 
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2 A Polynomial-time Heuristic algorithm 
2.1 The heuristic algorithm 

Let Rjj" be the set of non-negative real numbers and I the set of non-negative integers. 

Definition 1 Multi-constrained path problem (MCP): Given a directed graph G(V,E), a source 
vertex s, a destination node t } two weight functions W\ : E -> Rq and w 2 : E -> Rq , two constants 
d E Rj and c 2 G Rj ; the problem, denoted as MCP(G, s, t } w u w 2 , c u c 2 ), is to find a path p from 
s to t such that w\(p) < ci and w 2 (p) < c 2 if such a path exists. 

A path p which satisfies w^p) < c x and w 2 (p) < c 2 is called a solution to MCP(G, s, t } w u w 2 , ci, c 2 ). 
We assume that both weight functions are additive — the weight of a path is equal to the summation 
of the weights of all edges on the path. 

Definition 2 wx-weight and w 2 -weight: For a path p = v 0 -» vi Vk, 

k k 
W!(p) = E Vi) and w 2 (p) = .E ^(v.-i, t>i) 

™i(p) is called the wi-weight and iu 2 (p) the w 2 -weight of the path p. 

MCP(<3, 5, t, itfi, w 2 , ci, c 2 ) is NP-complete [12]. We provide a polynomial-time heuristic solution 
to this problem. The algorithm contains two steps: 

1. Create a new weight function w 2 : E — ► J. 

w' 2 {u,v) = (1) 

c 2 

where a: is a given positive integer. We reduce the original problem MCP(G, s, t } w u w 2 , Ci, c 2 ) 
to a new, simpler problem MCP(G, 5, t, w u ™ 2) c 1} x). 1 

2. Solve MCP(G, 5, t, w h w 2 , ci, x) in polynomial time. 

The algorithms for Step 2 will be discussed in Section 2.2. We assume for the moment that a 
solution of MCP(G, s, J, Wi, w 2) ci, a;) can be found in polynomial time if there exists one. 

^ote that the value of a: is chosen by the algorithm. It does not depend on the input values of G } a,tfli,it>2,ci and 
C2 . This is the essential reason for Step 2 of the algorithm to be solvable in polynomial time. As we will see shortly, 
a larger x means a higher chance to find a satisfactory path and a higher overhead. 
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Since MCP(G, s, t y w 1} w 2i ci, c 2 ) is NP-complete, we are not trying to find a solution for it when- 
ever there exists one. The idea is to reduce it to a simpler and solvable problem, MCP(G, s, t, wi, w 2l C\ t z), 
which has a "coarser resolution" — a? is a given finite integer and the range of w 2 is I. Theorem 1 
guarantees that a solution to the simpler problem must be a solution to the original problem. 

Theorem 1 A solution to MCP(G, s, t } wi, w 2l Ci, x) must also be a solution to MCP(G, s, J, iui, w 2 , ci, c 2 ). 

Proof: Let p be a solution to MCP(G, s, t, wi, w 2 , Ci, sc). Hence, < ci and w 2 (p) < x. In order 

to prove p is also a solution to MCP(G, 5, t } wi, w 2 , c x , c 2 ), it suffices to prove w 2 (p) < c 2 . By the 
definition of w 2} we have 

, r w 2 {u,v)x~" 

w ' 2 (u t v)= — s — by(l) 

c 2 

, w 2 (u,v)-x 

w' 2 (u,v)-c 2 
w 2 {u,v) < — 



Therefore, we have 



w 2(p)= w 2( u > v ) 

(u,v) on p 
(u,v) on p 

(u,v) on p 

< — • SB 



= c 2 



w 2 (p) < c 2 and hence the theorem holds. □ 

Corollary 1 Let P be the set of solutions to MCP(G, s, £, ti/i, tu 2 , Ci, c 2 ) and P' be the set of 
solutions to MCP(G, s } t y wi } w 2l ci, x). Then, 

P' CP 
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( a ) MCP(G, s, t, Wj , w 2 , 8.0, 20.0) ( b ) MCP(G, s, t, w, , , 8.0, 5) ( c ) MCP(G, s, t, w, , w^, 8.0, 10) 



Figure 1: (a) The path shown by the dotted line, 5— ► u -» v — > t, is a solu- 
tion to MCP(G } s,t 1 wi,W2 ) 8.0, 20.0). (b) Consider x = 5. The problem is reduced to 
MCP(s, t } tiii, w 2} 8.0, 5), which does not have a solution, (c) If x is increased to 10, the prob- 
lem is reduced to MCP(s, J, Wi t w 2l 8.0, 10), which has a solution, s — > u — > v — > t. Note that w 2 in 
(c) is different from that in (b) because the values of x are different. 

Corollary 2 Let P 1 be the set of solutions to MCP(G, s, t, w u w 2l c u x). The heuristic algorithm 
succeeds in finding a solution to MCP(G, s, t, wi, w 2l ci, c 2 ) if and only if P 1 / 0. 

The converse of Theorem 1 is not necessarily true — a solution to MCP(G, s, t, Wi, W2, ci, c 2 ) may 
not be a solution to MCP(G, s, t, W\ y w 27 ci, a:). Figure 1 gives an example. The original problem 
MCP(G, wi, t/>2j 8.0, 20.0) has a solution s — )> u — > v ->> t (Figure 1 (a)) . Suppose x = 5 and 
the problem is reduced to MCP(G, s, w\ } w 2} 8.0, 5) (Figure 1 (b)). The path s->u— »v-»£is 
not a solution to the new problem. In fact, there is no solution to the new problem. 

Hence, our heuristic algorithm may not find a solution for MCP(G, s, t } W2, c\ } C2) even when 
such a solution exists, because the solution set P 9 of the new problem MCP(G, 5, £, w u w 2l ci, x) can 
be empty. Fortunately, whether P r is empty or not is to some extend predictable and adjustable 
— by assigning a larger z, we have a better chance for P' to be non-empty. 

Theorem 2 2 Let a path p be a solution to MCP(G, s, t } Wi, w 2 , ci, c 2 ) and I is the length of p. If 

w 2 {p) < (1-— ) c 2 (2) 
x 

then p is also a solution to MCP(G, s, t, w\, w 2 , ci, a:). 

2 Theorem 2 can be rewritten as: A path p from a to t is a solution to MCP(G, s, t, w\ , tpj, ci, a;) if 

™x(p) < ci A ™ 2 (p) < (1 - ^) *c 2 
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Proof: Since p is a solution to MCP(G, s, t, w 2} Ci, c 2 ), we already have wi(p) < c\. In order to 
prove p is a solution to MCP(G, s, t y w u w 2i ci, 3), we only need to prove u> 2 (p) ^ 



^2(P)= £ ^K^) 
(«,v) on p 

= V r »»(".t>)-« n by(1) 

(u,v) on p 

< E (^^ + 1) 

(u,«) on p 

= f"" 13 U '2(«,U)+ J] 1 
(u,v) on p on P 

X 

= — -w 2 (p) + l 

C2 

<f (l-^) c 2 + Z by(2) 

C2 X 

= (x-l + l) + l 
= x + l 

Because both w 2 (p) and x are integers, w' 2 (p) < x. Therefore, the theorem holds. □ 

Theorem 2 means that if there exists a path p which is overqualified — not just w 2 (p) < c 2 but 
W2(p) < (1 — ^jp) * C2 — then after we reduce the original problem to MCP(G, 5, t } w\ } w' 2l ci, x), the 
new problem still has solutions (p is one of them). Hence, we can solve MCP(<7, 5, i, wi, it;^, ci, a;) 
to find a solution, which must also be a solution to the original problem as stated by Theorem 1. 

Corollary 3 Let P be the set of solutions to MCP(G, 5, t } u/i, w 2 , ci, c 2 ) and P' be the set of 
solutions to MCP(G, s, £, wi, ^3, Ci, 2). Then, 

P' / 0 if P ^ 0 and 3p G P,w 2 (p) < (1 - ^— ^) *c 2 
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Theorem 3 3 Let P be the set of solutions to MCP(G, s, t, w h w 2 , ci, c 2 ). The heuristic algorithm 
succeeds in finding a solution to MCP(G, 5, t, w u w 2l ci, c 2 ) if 

P / 0 and 3p e P, w 2 (p) < (1 - — ) • c 2 

x 

Proof: By combining Corollaries 2 and 3. □ 

Note that, by Theorem 3, if P ^ 0 and 3p G P, w 2 (p) < (1 — lj jr)c 2l then our heuristic algorithm 
is guaranteed to find a solution. However, even when P / 0 and Vp G P, w 2 {p) > (1 - ^)c2, the 
algorithm still has a good chance to find one, though that is not guaranteed. Figure 1 gives 
an example. In Figure 1 (a), MCP(G, s, w lf w 2 , 8.0, 20.0) has only one solution, which is p = 
5 -> w v -> £. w 2 (p) = 19.5. Consider x = 10 in Figure 1 (c). The problem is reduced 
to MCP(G,s,*,u;i,t4,8.0,10). Because i = 3 and c 2 = 20.0, (1 - ^)c 2 = 18.0. Hence, the 
condition w 2 (p) < (1 — ^)c 2 is not satisfied. However, p is still a solution to the new problem 
MCP(G, 5, t y wi t w 2) 8.0, 10). We thus conclude that w 2 < (1 — ^")c 2 is a sufficient but not a 
necessary condition for a solution of the original problem to remain as a solution of the new 
problem. 

How restrictive is the condition w 2 (p) < (1 — ^)c 2 ? Suppose a distributed multimedia applica- 
tion wants a connection in a network with delay no more than 300ms and delay jitter no more than 
100ms. Suppose there exists a path with a length of 5, a delay bound of 300ms and a delay-jitter 
bound of 90ms. Then with x = 40 our algorithm will definitely find a path for the application. 

Let us consider the worst case where the longest loop-free path in G is |V|. Table 1 shows the 
conditions under various different values of x. A larger x corresponds to a more relaxed condition, 
which implies a better chance to find a solution and a higher overhead — as we shall see in Section 
2.2, the overhead is proportional to x when an extended Bellman-Ford algorithm is used to imple- 
ment Step 2 of the heuristic algorithm. Take the case x — 10|V| as an example. The corresponding 
condition is w 2 (p) < 0.9c 2 . It means that, given an arbitrary problem MCP(G, s } t, w\, w 2) ci, c 2 ), 
if the problem has a solution p such that w 2 (p) < 0.9c 2 , then our heuristic algorithm is guaranteed 
to find a solution for the problem, provided a; is as large as 10|V|. 

3 Theorem 3 can be rewritten as: The heuristic algorithm succeeds in finding a solution to 
MCP(G, a, fc, ttfi , u>2>ci , ca) if there exists a path p from s to t such that 

wi{p)<d A w 2 (p) < (1 - -c 2 
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X 


w 2 {p) < (1 - ^) • c 2 


2\V\ 


w 2OO < 0.5c 2 


5\V\ 


w 2 (p) < 0.8c 2 


10\V\ 


*" 2(p) < 0.9c 2 


100|F| 


u; 2 (p) < 0.99c 2 


+00 


w 2 (p) < c 2 



Table 1: x v.s. tu 2 < (1 - ^) • c 2 



For real applications, our algorithm is expected to perform much better than what stated in 
Table 1. Consider the QoS routing problem by which our heuristic algorithm is originally motivated. 
The lengths of the routing paths between two nodes in the network should be much less than \V\. 
Therefore, to achieve a condition such as w 2 {p) < 0.9c 2 , the value of x can be far less than 10|V|, 
which means the actual overhead of the algorithm can be much lower than that of the worst case. 
Since the shortest paths are often preferred to save resources, a practical value for x can be 10d Stt , 
where d 8yt is the distance between 5 and t. 

One important question still remains: Is the new problem in Step 2 solvable in polynomial time? 
We answer this question in Section 2.2. 

2.2 The extended Dijkstra's and Bellman-Ford algorithms 

The problem MCP(G, s, t, wi, w' 2) Ci, sc), where w f 2 : E 6 I, is solvable in polynomial time for a given 
x G I. Note that no restriction has been placed on G, s, w\ and Ci, which can take any values 
allowed by Definition 1. 

An extended Dijkstra's shortest path algorithm (EDSP) and an extended Bellman-Ford algo- 
rithm (EBF) are presented in Figure 2 to solve MCP(G, s, t } Wi, w f 2l Ci, ar). An algorithm similar 
to EBF, in its distributed implementation, has been proposed by Jaffe [4]. We will discuss the 
difference between our algorithm and the Jaffe's algorithm in Section 3. 

For each vertex v G V and each integer k G [0..ar], a variable d[v, k] is maintained, which is an 
estimation of the smallest wi-weight of those paths from s to v whose tu^-weights are k. Let 

Min L {tui(p)} 

where P(v 1 k) = {p \ w f 2 (p) = A, p is a path from s to v}. The value of d[i/, fc], initially +00, is 
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Initialize(G, s) 




begin 


(1) 


for each vertex v E V[G], each i £ [Q..x] do 


(2) 


d[v, i] := oo 


(3) 


ir[v, i] := NIL 


(4) 


for each i £ [0..e] do 


(5) 


d[s,i]:=0 




end 




Relax(u, A;, v) 




begin 


(6) 


A;' := A; -f- w^u, v) 


(7) 


if Aj' < x then 


(8) 


if k'] > d[u, k] + w x (u, v) then 


(9) 


A;'] := d[uj k] + v) 


(101 






end 








Vk ACT 1 T1 
U 111 


fill 


Tn itiali "Z**! #■» c 1 
lili biailBCl Ur y a l 


(12) 


5 :- 0 


(13) 


Q:= {(u } k)\u£ V[G],ke [0..e]} 


(14) 


while Q ^ 0 do 


(15) 


find (u k\ G Q such that dta Jfel — Min -Tdki' jfe'H 




(16) 


Q:=Q-{(u,k)} 


(17) 


S:=S + {(u,k)} 




/* Note that the for loop iterates on different u's with fixed values for u and k. 




The values of u and Aj are changed by the outer while loop. */ 


(18) 


for each outgoing edge of u } (u, v) 6 E do 


(19) 


Relax(u, Aj, v) 




end 




EBF(G, 5) 




begin 


(20) 


Initialize(G, s) 


(21) 


for i := 1 to |F[G]|- 1 do 


(22) 


for A; := 0 to x do 


(23) 


for each edge (u, u) do 


(24) 


Relax(u, A;, v) 




end 



Figure 2: The Extended Dijkstra's Shortest Path algorithm (EDSP) and the Extended Bellman- 
Ford algorithm (EBF) 
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always greater than or equal to £(v,fc). During the execution, EDSP (EBF) makes better and 
better estimation and d[v } k] becomes closer and closer to, and eventually reach, S(v } k). 

When EDSP (EBF) completes, d[v, k] = S(v } k),v E V,k e [0..x]. There exists a solution, i.e. 
a path p from s to t such that w\{jp) < C\ and w f 2 (p) < x, iff 3k G [0..z], d[t, k] < C\. The path is 
stored by 7r. The variable 7r[J, k] keeps the immediate preceding vertex (called predecessor) of £ on 
the path. Hence, the path can be recovered by tracing the variable ir of all intermediate vertices 
till reaching the source s. 

Two additional variables, S and Q, are required by EDSP. 

5 = {(v, Jb) | d[v, k] = 6{v, Jfe), v e V, k G [0..*]} 

Q = {(v, k) | d[v, k] > 6(v, k), v G V, k e [0..as]} 

where the notation (v, fc) simply means a pair of values, v G V and A: G [0..x]. Initially, 5 = 0 and 
Q = {(u, A)|t; G V, Aj G [0..a:]}. By the while loop (lines 14-19 in Figure 2), each iteration moves a 
pair from Q to S and adjusts the wi-weight estimation by calling Relax(u, k, v). When Q — 0, the 
algorithm completes. 

A more detailed presentation of the original Dijkstra's and Bellman-Ford algorithms, which our 
algorithms are based on, can be found in [2]. 

The time complexity of ESDP is (x 2 V 2 ). The maximum size of Q is (x + 1)V. Hence, line 15 
can be done within 0(xV) . There can be at most (x + 1)V iterations of the while loop and thus 
the total time for line 15 is 0(x 2 V 2 ). The for loop of lines 18-19 has {x + 1)E iterations in total 4 
because Relax(u, A:, v) is called once for every (u, x;) G E, k G [0..x]. In each iteration, Relax(u, k } v) 
takes O(l). Hence, the time complexity for this part is 0(xE). The total time complexity is 
0(x 2 V 2 + xE) = 0{x 2 V 2 ). The time complexity of EBF is 0(xVE), because line 23 is executed 
for at most (x + 1)(V — 1)E times. The space complexities of both algorithms are 0(xV). 

Let us consider the time complexity of our heuristic algorithm in Section 2.1. Step 1 of the 
algorithm takes 0(E). Step 2 of the algorithm is implemented by EDSP or EBF. Therefore, the 
total time complexity is 0(x 2 V 2 ) when EDSP is used or 0(xVE) when EBF is used. The time 
complexity is polynomial because the value of x is given by the algorithm. For example, if we let 
x = 10\V\ and use EBF in Step 2, the time complexity is 0(V 2 E). 

4 (z + 1)E iterations are the combination result of the outer while loop and the inner for loop. The while loop 
iterates on u and fc, and the for loop iterates on v. 
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2.3 Algorithms for more than two constraints 

We have been studying MCP with two constraints so far. In this section, we shall generalize the 
algorithms for more than two constraints. We define MCP(G, s, t, W\ % w n) c x , c n ), n > 2, as to 
find a path p from s to t such that W{ < Cj, for each i < n. 
The heuristic algorithm is generalized to be: 

1. For each i E [2..n], create a new weight function : E — > I. 

r w i (u i v)x i ~ } 
Wi(u,v)= 

where X{ is a given positive integer. We reduce the original problem MCP(G, s, t y w\ } w ni c\ , ... 
to a new, simpler problem MCP(G, s.t.w^w^ c x , a? 2 , ^n), which is to find a path p 

such that wi(p) < ci and w f { (p) < X{, 2 < i < n. 

2. Solve MCP(G, s, t, w u w' 2 , w' n) c u x 2 , x n ) by EDSP or EBF. 
The generalized EDSP and EBF can be found in Figure 3. 

The time complexities of EDSP and EBF are 0{x\ • ... *x\ V 2 ) and 0(x x '...'X n VE) } respectively. 
For sparse graphs in which E = 0(V) } EBF is likely to perform better than EDSP, especially when 
n is large. The time complexity of the generalized heuristic algorithm equals that of EDSP or EBF, 
depending on which is used for Step 2. 

2.4 Multi-Constrained Routing 

Multi-Constrained routing is an important application of MCP. A network can be modeled as a 
graph G(Vj E). V is the set of routing nodes and E is the set of communication links. Each link has 
two properties: delay and cost. 5 The delay (cost) of a path is the summation of the delays (costs) 
of all links on the path. Given a source node s and a destination node £, the multi-constrained 
routing problem is to find a path p from s to t such that delay(p) < D and cost(p) < C, where 
D and C are the required end-to-end delay bound and cost bound, respectively. This is clearly a 
MCP problem. The routing algorithm is presented below. 

5 The cost of an edge can be measured in dollars, or it can be a function of a given system metric such as bandwidth 
utilization or buffer utilization. 
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Initialize(G, s) 
begin 

for each vertex v G V[G], each k{ G [CzJ, 2 < % < n do 

d[v f A 2 f —i *n] : = co 

7r[v,* 2 ,— i An] := i^/Z 
for each At £ [0..a?t]» 2 < i < n do 
d[s, A 2 , A n ] := 0 

end 



Relax(u, A 2 , A n ,v) 
begin 

for each i G [2..n] do 

A< := ki + w'ifav) 
if A^ < Zi, for each i G [2..n] then 

if A 2 , A n ] > d[u } A 2 , A n ] + v) then 
d[v, A 2 , AJJ := i[u, A 2 , A n ] + to^u, v) 
7r[v, A 2j AJJ := u 

end 



EDSP(G, 5) 
begin 

Initialize(G, 5) 
5:=0 

Q ~ {(u, A 2 , A n )|u G V[G], Ai € [0..**], 2 < i < n} 
while do 

find A 2 , A n ) G Q such that <i[u, A 2 , A n ] = Mm {d[t*', A 2 , k' n ]} 

G:=0 -{<!*, A 2 ,...,A n )} 
S:=S + {(u } A 2 ,...,A n » 
for each edge (u, v) G ^ do 
Relax(u, A 2 , A n , v) 

end 



EBF(G, 5) 
begin 

Initialize(G, s) 
fori:=lto |V[G]|-ldo 

for each A< G [0, asj, 2 < i < n do 
for each edge (u, do 
Relax(u, A 2 , A n , v) 

end 



Figure 3: The generalized EDSP and EBF 
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1. Create two new functions new-delay : E -> I and new-cost : E -» J. 

new-delay[u, v) = — 

/ x r costlu, v) • x n 
neiu-cosi(u, v) = — 

where x = coef X d, it , coe/ is a given positive integer and d B)t is the distance from s to 
t. We reduce the original problem MCP(G, s,t, delay, cost, D,C) to two simpler problems, 
MCP(G, s, t, delay, new-cost, D, x) and MCP(G, s, t, new-delay, cost, x, C). 

2. First, solve MCP(G, s, t, delay, new-cost, D, x) by EDSP or EBF. If a solution is found, return 
the found path and terminate; otherwise, solve MCP(G, s, t, new-delay, cost, x, C). 

A distributed implementation of the above algorithm is possible. However, without further 
simplification of the algorithm, the number of messages to be sent is expected to be very high 
[4]. Hence, we assume a source routing strategy, which was also adopted by routing algorithms in 
[6, 5, 8, 12]. It requires every node to maintain the state information of the network, which can be 
done by the link-state algorithm [7]. The routing path is determined locally at the source node. 

The proposed routing algorithm applies the heuristic algorithm (Section 2.1) twice, reducing 
delay and cost to new-delay and new-cost, respectively. Hence, it guarantees to find a solution 
when either of the following two conditions is satisfied by a path p from s to t (see Theorem 3): 

1. delay(p) < (1 — ^) • D A cost(j>) < C (Heuristic condition one) 

2. delay(p) < D A cost{p) < (1 — ^) • C. (Heuristic condition two) 

The above two conditions are called the heuristic conditions. Note that we use the heuristic 
conditions to theoretically study the performance of the algorithm but they are not part of the 
algorithm. 

The heuristic conditions are sufficient conditions in the sense that if one of them is satisfied 
our proposed routing algorithm guarantees to find a solution — it makes the performance of the 
algorithm predictable to some extend. However, as we discussed in Section 2.1 and showed by 
the example in Figure 1 (c) , they are not necessary conditions — even when none of the heuristic 
conditions is satisfied, our algorithm may still find a solution. Hence, the heuristic conditions give 
a pessimistic estimation (or lower bound) of the performance of the proposed routing algorithm. 
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Figure 4: the topology of the experimental network 



2.5 Experiments 

We know from the routing algorithm proposed in Section 2.4 that x = coef X d Bjt . What is the 
relationship between coef and the performance of the algorithm and how large should coef be? 
We answer the questions by simulation. 

The network topology used in our simulation is shown in Figure 4, which expends the major 
circuits in ANSNET [1] by inserting additional links to increase the connectivity. For each routing 
request, the values of s } t, delay, cost, D and C are randomly generated. The delay values of 
the links are uniformly distributed in the range of [0..50ms], and the cost values of the links are 
uniformly distributed in [0..200]. The performance metric we considered was success ratio. 

number of requests successfully routed 

success ratio = : 

total number of routing requests 

We studied the success ratio with respect to coef, D and C. The larger the value of coef, the 
higher the probability for the heuristic conditions to be satisfied, which leads to a higher success 
ratio. The smaller the values of D and C, the tighter the constraints of a routing request, which 
leads to a lower success ratio. 

The experiment results are presented in Figures 5-9. The x axis represents coef and the y axis 
represents the success ratio. The dimensions of D and C are shown by different figures. Let us 
take Figure 5 as an example. Each point in the figure is taken by running one thousand randomly- 
generated routing requests. The values of D and C of all the requests are uniformly distributed 
in [100,115ms] and [400,460], respectively. For the purpose of comparison, we implemented an 
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The average source-destination distance of all connection requests is 3.2 hops. 
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the optimum algorithm 
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percentage of requests which satisfies one heuristic condition b— 
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23456789 10 
coef (Note that the value of x is equal to coef times the source-destination distance.) 

Figure 5: Success ratio with respect to coef when D £ [100, 115ms] and C £ [400,460]. Note that 
x = coef • d 3)t . The average distance from the source to the destination of all experimental routing 
requests is 3.2 hops. 

optimum algorithm, which searches all possible paths for a solution with an exponential time 
complexity. There are three lines in the figure. The highest horizontal line shows the success ratios 
of the optimum algorithm. The lowest line shows the percentage of requests which satisfy one of 
the heuristic conditions, 6 Note that the lowest line is calculated purely based on the heuristic 
conditions. It gives a lower bound of the performance of the proposed routing algorithm. The 
actual running performance of the algorithm is given by the middle line, which is constantly above 
the lowest line. It means that the statistical performance of our algorithm is better than what 
predicted by the heuristic conditions. Figures 6-9 present the success ratios when D £ [75, 90ms] 
and C £ [300,360], D £ [125,140m*] and C £ [500,560], D £ [150,165ms] and C £ [600,660], 
D £ [150, 165ms] and C £ [600, 660], respectively. Larger values for D and C result in more relaxed 
delay and cost constraints and thus higher success ratios as seen from the figures. 

The simulation results show that the success ratio of the proposed routing algorithm approaches 
that of the optimum algorithm when coef is increased. With coef > 4, the performance of our 
algorithm is close to that of the optimum algorithm. 

6 The heuristic conditions provide a theoretically-provable performance estimation for the algorithm (see Theorem 
3). They are useful as a guide line in the design of a network which requires a guaranteed routing performance. 
However, please be aware that they are not part of the routing algorithm. 
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Figure 6: D G [50, 65ms], C G [200, 260] 



Figure 7: D G [75, 90ms], C G [300,360] 



Tho average source-destination distance of all connection requests is 3.2 hops. 



The average source-destination distance of all connection requests is 3 2 hops. 
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Figure 8: D G [125, 140ms], C G [500, 560] Figure 9: D G [150, 165ms], C G [600, 660] 
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3 Related Work 



Much work has been done in QoS routing recently [5, 9, 10, 11, 12]. Some routing algorithms 
consider a single constraint. Plotkin discussed the competitive routing strategy in [9], which con- 
siders only the bandwidth requirement. Salama et al [10] proposed a distributed algorithm for the 
delay-constrained routing. The Sun's algorithm [11] improved the worst-case complexity of the 
Salama's algorithm. Though both algorithms use heuristic approaches trying to minimize the cost 
of the found route, the cost is not required to be bounded. 

The multi-constrained routing was studied in [5, 12]. Wang and Crowcroft [12] uses the Dijk- 
stra's algorithm in their bandwidth-delay-constrained routing, i.e., finding a path whose bandwidth 
is not less than a given lower bound and whose delay is not greater than a given upper bound. 
Note that the bandwidth-delay-constrained routing is not a MCP problem by Definition 1. Both 
weight functions in MCP, Wi and w 2 , are additive — the wi-weight (w 2 -weight) of a path is 
equal to the summation of the wi-weight (w2~weight) of all edges on the path. In the bandwidth- 
delay-constrained routing, however, one of the weight functions, bandwidth, is not additive — the 
bandwidth of a path is equal to the minimum bandwidth of all edges on the path, which makes 
the problem much easier and can be solved in polynomial time. 

Ma and Steenkiste [5] showed that when a class of WFQ-like (Weighted Fair Queuing) schedul- 
ing algorithms are used, the problem of finding a path satisfying bandwidth, delay, delay-jitter, 
and/or buffer space constraints is solvable by a modified version of the Bellman-Ford algorithm 
in polynomial time. The reason is that when the WFQ-like scheduling algorithms are used, the 
metrics of delay, delay-jitter and buffer space are no longer independent from each other and all 
of them become functions of bandwidth, which simplifies the problem and makes it solvable in 
polynomial time. 

All the above algorithms can not solve MCP whose weight functions are assumed to be inde- 
pendent. The work closest to ours was done by Jaffe [4]. Jaffe proposed a distributed algorithm 
solving MCP with a time complexity of 0(V 5 b log Vb) y 7 where b is the largest weight of all edges 
in the network. The complexity is pseudo-polynomial because the run time is polynomial in 6, the 
largest number in the input. See [3] for the definition of the pseudo-polynomial time complexity. 

A heuristic algorithm with a polynomial time complexity was also proposed in [4] to approximate 

7 The Jaffe's algorithm finds a solution for every pair of nodes in the network whereas our heuristic algorithm does 
that for a single pair. 
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Figure 10: The two constraints are Wx(p) < 8.0 and w 2 (p) < 8.0. The path, s — > u — ► t, minimizes 
wi(p) + W2{p) but is not a solution. The path, s — > v — > t, does not minimizes wi(p) + W2(p) but 
is a solution. 

the MCP problem. The algorithm, referred to as the Jaffe's approximation algorithm, finds a path 
which minimizes wi(p) + d • w 2 (p), 8 where d is a given constant. Clearly, it can be implemented by 
either the Dijkstra's or the Bellman-Ford shortest path algorithm. However, minimizing wi(j>) + 
dw 2 {p) may not lead to a solution of MCP. Figure 10 gives an example for the case 6=1. The 
problem is MCP(G, s, t } ti/i, w 2i 8.0, 8.0). The path, s — > u J, minimizes wi(j>) + w 2 (p) but is 
not a solution because the constraint w 2 (p) < 8.0 is not satisfied. The path, s — ► v — > t, does 
not minimize wi(p) + w 2 {p) but is a solution. Hence, the JaflFe's approximation algorithm fails in 
finding a solution in this example. However, our routing algorithm proposed in Section 2.4 finds 
the solution with x = 2. 

The time complexity of the Jaffe's approximation algorithm is 0(V 2 ) when the Dijkstra's al- 
gorithm is used or 0(VE) when the Bellman-Ford algorithm is used. It is much better than the 
complexity of our heuristic algorithm, which is 0(x 2 V 2 ) when EDSP is used or 0(xVE) when EBF 
is used. Therefore, in practice, we can apply the Jaffe's approximation algorithm first. If it fails in 
finding a solution for MCP, our routing algorithm is executed next. 

4 Conclusion 

Any multi-constrained routing problem which involves two or more additive weight functions such 
as delay and cost is NP-complete. We formalized it as an MCP problem (Definition 1) and proposed 
a heuristic algorithm with a polynomial time complexity. The algorithm first reduces the problem 
MCP(G, 0,t,tui,t!/2,ci,C2) to a simpler one MCP(G, s, t, W\ } w 2 , Ci, x) } and then uses an extended 
Dijkstra's (or Bellman-Ford) algorithm to find a solution for the new problem in polynomial time. 

8 The paper [4] used different terminologies and notations from ours. wi(p) + d-W2(p) was written as L(p) + d-W(p), 
where L and W correspond to w\ and u>2i respectively. 
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We showed the correctness of the algorithm by proving that any solution found for the simpler 
problem must also be a solution to the original problem. We showed the effectiveness of the 
algorithm by proving that the simpler problem must have a solution if the original problem has 
a solution p and w 2 {jp) < (1 — ^) c 2 ) where I is the length of p and x is an integer given by 
the algorithm. With an increasing x, the condition w 2 (p) < (1 - ^) c 2 is gradually relaxed and 
approaching the original constraint, w 2 (j>) < c 2 . Even when the condition w 2 (j>) < (1 — lj ^-)c 2 is not 
satisfied, our algorithm still has a good chance to find a solution because w 2 (p) < (1 — ^)c 2 is a 
sufficient but not a necessary condition. The statistical performance of the heuristic algorithm was 
studied by experiments, which showed that higher performance of the algorithm can be achieved 
at the expense of higher overhead. 
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